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I. INTRODUCTION 


In this paper we consider the production of a pair qq of heavy quarks close to threshold; 
that is to say, for small values of |u| where 

v = \J 1 — 4 m 2 / s . (1.1) 

Here m is the quark mass, and s 1 / 2 the center of mass energy of the qq pair. Above threshold 
v coincides with the velocity of either quark; but we also study the production of qq bound 
states. 

We take the quarks to be heavy, so that we may be able to apply a perturbative QCD 
analysis to them. Thus we study production of cc, bb and tt. (For the last case, and in this 
first paper, we will neglect the effects of t decay). With a view to applications to production 
by e + e~ collisions, we take the qq to be produced by either a vector V or axial A current^;: 

V(x) = q(x) q(x) 

A(x) = q(x) 7^75 q(x). ( 1 . 2 ) 

Sum over omitted color indices is understood. We will then study the correlators, 

n \ iU {p) = if d 4 x e ip ' x (vac\ TV^ix) 14(0) | vac) 

= ( ~P 2 gvv + P^Pu) n v(p 2 ) , (1-3) 

n f iU {p) = i J d 4 x e w ' x (vac\ T A^(x)A u ( 0 ) | vac) 

= ( -p 2 g^ + p^Pu) n A {p 2 ) + p»Pv n P (p 2 ), (i.4) 

where | vac) denotes the physical vacuum. As is known, the production cross section for 
e + e _ —> qq may be straightforwardly written in terms of the imaginary parts of the correla¬ 
tion functions 

imn(p 2 ) , n = n V:AP ; (1.5) 

1 Our results may be extended with little effort to scalar or pseudoscalar correlators. 
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only Ily, HU give sizeable contributions, and they are the quantities we will study here. 

To lowest order in perturbation theory (he., the parton model) we may neglect the 
interactions of qq. The II are then obtained with a simple one-loop evaluation and one has 


Im4 0) (s) 

Imhfts) 


N c v(3 — v 2 ) 
127T 2 


6{s- 


N c 
12vr 


v 3 9(s 


4 m 2 ) 



(1.6.a) 
(1.6.b) 


Here N c = # colors = 3, and the subscript zero in the n indicates that the strong inter¬ 
actions are neglected. Of course strong interactions are most important near threshold: its 
incorporation is precisely the subject of the present paper. 

This article is organized as follows. In section 2 we give the expression for the contri¬ 
bution of the qq bound states to the ImE This, we hope, will serve to clear some of the 
misunderstandings found in the standard literature. We also deal, in this section, with the 
Imll(s) above threshold, but in the nonrelativistic regime (u 2 <C 1), explaining the use of 
the final state interaction theorem to incorporate strong interactions, which again should 
clarify some of the existing fog. 

In section 3, we use the known evaluations of wave function for bound states, and the one 
obtained here in the continuum, to give explicit formulas for ImII(p 2 ) below threshold and 
above. In this last case radiative and nonperturbative contributions are included for the first 
time. The article is concluded in section 4, where we discuss in detail the important case of 
ImHy(p 2 ) above threshold and the contribution of this region to the evaluation of Hy(g 2 ) 
for q 2 < m 2 . An Appendix is also provided for the technical details of the calculations. 


II. THE IMAGINARY PART OF II AROUND THRESHOLD 

A. Bound state contributions to Im II. 

We will carry over the detailed calculations for Im n^; then we will indicate the corre¬ 
sponding results for Im n^. If we have a bound state of qq with momentum k, k 2 = M 2 , 
and third component of spin A, normalized to 
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(k, X | k',X') = 2k Q 5 xx ,5(k-k') , 
then its contribution to Im U^fip) is 

Im n|2’”' e (p) = i J d‘x e ‘ P ’ I Xi / JjT (vac\J^(x) |fc,A) (k, A| J„( 0 ) |t>ac) 

= i(27r)“5(p 2 - M 2 ) ^ V„(p, A) V„*(p, A) . (2.1) 

Z A 

In the above equation, 

Vn(p, A) = (uac| g( 0 ) 7 ^g( 0 ) \p, X) (2.2) 

is the (p-space) wave function of the bound state by definition. The connection with the 
rc-space wave function may be carried over immediately in the nonrelativistic limit. In the 
c.m. referencial one thus finds, for e.g. X = +1, and with k = pi + p 2 , 

(vac\ g( 0 ) 7 M g( 0 ) \k, X = 1) 

v(pi, 1/2) 7m u(p-2, 1/2) T(0) . (2.3) 

Here T(0) is the x-space wave function evaluated at r = 0: thus, only states with £ = 0 
contribute. 

It is perhaps not idle to note that Eqs. (2.1, 2.3) are exact, the last if taken as a definition 
of 'I'(O). Substituing one into the other, we finally obtain 

Imn p v ole (p 2 ) = ^ 5{p 2 - M 2 ) |i?o(0)| 2 , (2-4.a) 

and R 0 = (4 7 t) 1//2 T f=0 is the radial wave function. 

For the axial correlator the evaluation is slightly more complicated because the orbital 
angular momentum of the bound states is t = 1. One finds, 

im n yv) = hp 2 - m 2 ) |«;<o)| 2 , 

R'fir) = dRi(r) / dr ; (2.4.b) 


VK I k 0 

(2 7 r ) 3 / 2 \l 2p w p 20 


Ri is the radial part of the £ = 1 wave function. In the nonrelativistic limit the Rg are 
normalized to 
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B. Imll above threshold. 


We will give a detailed discussion for the case where we have the vector correlator. More¬ 
over, we will consider only the production via a virtual photon, neglecting the contribution 
of the Z. (This only for ease of discussion; the results will be valid quite generally.) In these 
circumstances the quantity Im H^ LV may be considered, for p 2 > 4 m 2 , to be proportional to 
the square of the production amplitude 7 * —► qq : 

Im ~ I (qq\ S ( 7 *) | 2 . 

We will consider that the interactions involved in this process are two: the electromagnetic 

Hie m = eQ f j d 3 x q(x)^^q(x) A fl (x) , 

and the QCD interaction described by a Hamiltonian Hjqcd that will be specified later. The 
final state interaction theorem then asserts that the amplitude (qq\ S |y*) may be evaluated, 
to first order in Hi em , but to all orders in Hiqcd, by means of the expression 

(qq\S\r) = i(*| , (2.5) 

where |T) is an exact solution of the Lippmann-Schwinger equation for qq states subject to 
the strong interaction: 

U) = U (0) ) + tt— rr Hiqcd I'f) . (2.6) 

— H o 

(In our case, and because only one wave contributes, we need not specify the boundary 
conditions in Eq. (2.6)). 

There is now a complication, as compared to the bound state case: |T) may now contain, 
besides qq , qq + n gluons. To the order we will be working and in the nonrelativistic 
approximation , the states | qq + n gluons) may however be neglected. The reason is that the 
amplitude for radiation of a gluon by heavy quarks is proportional to the velocities v q + v q , 
so that the contribution of these processes to Im will be of order (v q + ^-) 2 - ^ 
i.e, of the order of the relativistic corrections, which we are neglecting. Therefore, in the 
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nonrelativistic limit the state |\&) may be considered to consist only of qq, and can thus be 
represented by a wave function. We then obtain 

Im n v(s) = I R k o(0) I 2 Im n P(s) , (2.7.a) 

and also 

Im U A (s) = \3R' kl (0)/mv | 2 Im nj } (s) . (2.7.b) 

Here R k t is the radial part of the continuum wave function with k = mv, and normalized to 

too 9 7r 

J o dr r 2 R* ke (r) R k , e (r) = — S(k - k') . (2.8) 

The Im n (0 ^ are as given in Eq. (1.6). Again we would like to emphasize that Eqs. (2.7) are 
exact (in the nonrelativistic limit); approximations enter when evaluating the R k e, which 
will be the subject of the next section. 

III. THE qq WAVE FUNCTION CLOSE TO THRESHOLD 

A. The QCD interaction for heavy quarks at short distances. 

It has been known for a long time that the short distance interactions of a pair of heavy 
quarks may be described by perturbation theory; the leading nonperturbative corrections are 
then implemented taking into account the nonzero values of quark and gluon condensates 
in the physical vacuum: 

(qq) = (vac\ : q(0 ) g(0) : | vac) , 

(a s G 2 ) = a s (vac\ : G^ w (0) (7^(0) : | vac) . 

The details may be found in the classical SVZ papers^ where the correlators n /)jy are directly 
studied using the operator product expansion techniques; or the work of Leutwyler^ and 
Voloshin^ where the Green’s function method is employed to study the qq bound states. 
In particular, these last authors explicitly prove that no potential may describe the short 
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distance nonperturbative corrections to the qq spectrum and wave function; instead, one 
has to use an effective interaction which in the nonrelativistic limit is given by 

HiNP = —gr£ a t a , (3.1) 


with £ the chromoelectric field. One then takes, 

(vac | £ | vac) = 0 , (vac\ £ 2 \vac) ~ ( a s G 2 ) . 


(In our case the contribution of the quark condensate is negligible). From a practical point 
of view it has been made apparent in the detailed evaluations of Ref. [4] that a calculation of 
qq states, based on perturbation theory, and supplemented by the leading nonperturbative 
corrections, as given by Eq. (3.1), yields an excellent, essentially parameter-free description 
of the bound states of cc with n — 1 (n being the principal quantum number) and of bb 
states with n — 1,2 ; £ = 0,1. As was already known from the work of Refs. [2,3], the 
analysis breaks down for higher excited states where nonperturbative contributions get out 
of hand and calculation from first principles becomes impossible. This occurs for n > 1 in 
cc , and n > 2 for bb . For tt the distances involved are so short that a rigorous calculation 
becomes possible up to n ~ 5. 

Besides the nonperturbative contributions described by Eq. (3.1) we require also the 
interaction deduced using perturbation theory. In the nonrelativistic regime, and including 
one-loop correction^], we have the Hamiltonian 


H = H?} f + H 1 , 

u~(o) 1 A C F 5 S (/1 2 ) 

tL e ff — A 

jj m r 

„ c F fi a s lnr g 

11 1 = — - - ; 

27 T r 

here a s includes part of the radiative corrections, 


(3.2) 


2 The radiative corrections to the qq potential have been obtained by a number of authors; cf. Ref. 
[5] and Ref. [4] where they are completed and summarized. 
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a s (n 2 ) = a s (p 2 ) 


1 + ai + 1e ^ o/2 a s (/i 2 ) 


7 r 


_ IIC 4 - 4T f"/ , _ 310,-20^^ „_4 1 

« 3 ■ “1 36 ’ °r 3. k. 3 , J F 2 . 

and rif is the number of quark flavors with masses much smaller than m. The reason 
why we include part of the radiative corrections in is that this Hamiltonian, being 
Coulombic, may (and will) be solved exactly, whereas H l and H INP have to be incorporated 
in perturbation theory. 

A last point about Eq. (3.2) is the meaning of the parameter m there. As follows from 
the analysis of Ref. [4], this m has to be interpreted as the pole mass. That is to say, if S(tf) 
is the quark propagator, in perturbation theory, then m is such that A -1 (ft = m ) = 0. This 
m may be related to the MS mass, m(m 2 ) through the formula) 6 ! 


m = m ( m 2 ) \ 1 + C F a s( m ) ^ _ 2 C F ) 


n 


a s ym- 


n 


+ ... 


(3.3) 


Here K ~ 13.5 (an exact formula for K may be found in Ref. [6]). As for the numerical 
values of the masses, the analysis of Ref. [4] gives 


m c = 1570 ± 60 MeV ; m b = 4906 ± 85 MeV , 


(3.4) 


which corresponds to the MS masses 


m c (m 2 ) = 1306 ± 40 MeV , m b (ml) = 4397 ± 40 MeV , 


of course compatible with (but more precise than) the values obtained with the SVZ 
method! 1,7 ! 


m c (m 2 ) = 1270 ± 50 MeV , m b (m 2 ) = 4250 ± 100 MeV . 
We will thus take Eq. (3.4) as our input. For the top quark we elect to choose 


rn t = 165 ± 15 GeV 


( 3 . 5 ) 



B. The bound state wave functions. 


The energy levels calculated with Eqs. (3.2), (3.1) are^ 

C r®s(p 2 ) 2 C F (5 0 a% 


E ne = 2 m 


8 n 2 8 7T n 2 

Tte n e n 6 (a s G 2 ) ) 

2 (mC F a s ) 4 } 


In ■ 


fin 


mC F a s 


+ ip(n + £ + 1 ) 


(3,6) 


Here the e are numbers of order unity: 

_ 624 _ 1051 _ 769456 _ 9929 

610 “ 425 ’ 620 ““ "663” ’ 630 “ 463239 ’ 621 “ 9945 ’" ' 

other e n £ may be found in Ref. [4], and an analytic expression in Ref. [2], For the wave 
functions, the details are given in Ref. [3] and, especially, in the second paper of Ref. [4], 
One has 


Rntir) = 


n 2 a(n,£) 3 / 2 \J (n + £)\ 
Pnt = 2r/a{n,l) ; 


(n-£- 1 )! 


pie e M/2 (Pne) , 


(3.7) 


a(n, £) = 


mC F a s (n 2 ) 


In (nn/mCpOts) + i>(n + £ + 1) - 1 _ 

1 - 2 ?- 


This includes the radiative corrections, as obtained using Eq. (3.2). The full wave function 
at the origin, including nonperturbative contributions (Eq. (3.1)) is 


Rni{ 0 ) 


^1 + 5np(u, £)^j Rni{ 0 ) 


(3.8) 


and, for the first n, £, 


5np( 1 , 0 ) = | 
8np( 2 , 0 ) = | 
8np( 2 , 1 ) = | 


2 968 968 576} n(a s G 2 ) 

425 + 541 875j m 4 (C F a s ) e ’ 

3 828 736 753 025 024) n(a s G 2 ) 

1989 + 1318 707 J m 4 (C F a s f ’ 

3 299 840 33 026 904 064) n (a s G 2 ) 

1989 + 98 903 025 J m 4 (C F a s ) e ' 


Higher S^p(n, £) may be found in the second paper of Ref. [4], 
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C. The qq wave function in the continuum. 


The calculation of the R^e with the Hamiltonian of Eq. (3.2) is far from trivial. It is 
described in some detail in the Appendix. We present here the results for the modulus 
squared, at r — 0 : we have 


fl«(0) =[1 + 2 c,(k)] R k ,( 0) 


(3,9) 


Here R is evaluated with so that 


Rk o(0) 


nC F a s jv 

^ _ g-7T C f OIb/v ’ 


(3.10.a) 


3Rki(0)/mv 


= 1 + 


C F a s \ 7 r C F a s jv 


(3.10.b) 


Av 2 ) i _ e~^ c F ° s ! v 

a s = a s (/a 2 ) given after Eq. (3.2). The functions c^(k) are plotted in Figs. 1, 2 for i = 0,1, 
and displayed in detail in the Appendix. For t = 0,1, and small velocities, 

/% 


Co(fc) = 

ci(k) = 


47T 

A a s 

2 7T 


, iia 
In — 


( ka ) 2 ( ka ) 4 

1-2 7b + ^^ + E_E + ... 


12 


40 


, V a o 

In — - 2 7 E 


\ 4Q 

+ (ka) 2 - (ka) 4 J - 2 + — (ka) 


' 2 f 


2 7T 


[ka — (ka) 3 + ( ka ) 5 ) 


+ 


a 


in C F a s 


(3.11) 


A remarkable property of Eq. (3.11) is that, as k —» 0 the relevant scale for a s (/a 2 ) is 
not /a ~ k (as one would naively gness, and as assumed for instance in Refs. [8,9]), but 
H ~ 2 /a = mC F a s : the interaction saturates. 

The full wave function is obtained by adding the (leading) nonpertnrbative contributions. 
These may be deduced from the the SVZ calculations^. One finds 


\R k i(0)\ 2 = (l + 2S^ p ) \R ke (0) 


(3.12) 


For i = 0 one has 
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cNP _ ^ 

fc0 “ _ 192 m 4 v 6 


(3.13) 


Just like for the bound states, the nonperturbative correction blows up at threshold (v —> 0). 


Clearly, the calculation ceases to be valid when 
he., for a critical velocity v cr u such that 

( n 2 (a s G 2 ) 

Vcrit ~ V 192 A 


NP 


2\ \ VS 


is of the order of magnitude of |c^|, 


(3.14) 


IV. NUMERICAL RESULTS 


A. Production of qq above threshold. 


We now consider the quantities 


Rq(s ) = 12 7T Q 2 Im fly (s) , 


Rfis) = 12 7 T Q 2 Im n^(s) . 


(4.1.a) 


is essentially the ratio (e + e —> 7 * —» qq)/(e + e —> 7 * —> /x + /i ). To show clearly 
the various contributions, we will plot the zeroth and first order (in a s ) expressions for the 
R q (s), cf. Eq. ( 1 . 6 ): 


= LjV , 

RV , (1, W =„ N ‘Q 2 , {x “ V + IT + " j c ' a ‘ ( " 2 ) • 

R 'i (1>(s) „“i {I + i (1 “ + (i 1,1 ru; ~ |) - '"> 2 + ■ ■ ■} 


(4.2.a) 


C F a s (n 2 


n 


R^°\s) = N c Q 2 q v 3 , 

^ (1) ( S ) = |7tC f a s (fl 2 ) V 2 + . . .} . 


(4.2.b) 


The expression for R q ^ is actually known exactly^, but Eq. (4.2.a) is accurate enough 
for our purposes. Only the leading term (in v) in R^ ^ is known, and it is reported in 
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Eq. (4.2.b). It is worth noting that the Shwinger interpolation! 10 ! to i? 1 d), used by some 
authors, is not accurate enough for our purposes. In fact, it reproduces correctly the values of 
i? 1 d)(s) for v = 0, v — 1, but not the derivatives, i.e., the terms in v, v 2 ; (1 — v), (1 — v) 2 . 
For this reason we take the full (4.2.a) (see below, e.g. in (4.3.a)). 

Together with these perturbative evaluations we also plot the “exact” expressions ob¬ 
tained from Eqs. (3.9-3.13): 


n F, . Ar o f v(3 — V 2 ) ( 6 V 3 TT v 2 . „ 

R , w =„ N °<¥, + — + — I c - a - 


i - 


x [1 + 2 Co(k)] 


7T CpOLs/v 


n(a s G 2 ) ' 

192 m 4 v e 
\ / 

Oi s Q? s (/i^) , 


_ g-7T CpCts/v 

we have included the known 0(a s v, a s v 2 ) corrections, cf. Eq. (4.2.a). For R£ 

C 2 a 2 \ 7r C F a s /v 


(4.3.a) 


Rf(s) = N c Q 2 v A (1 + 2 Stf) [1 + 2 Cl {k)} 1 + 


F s 


1 /L 1 r 1 4v 2 I 1 _ e -FC F Os/v 


(4.3.b) 


The quantities 5^ may be found in Ref. [11]. The calculation of R 1 will be made for 
q — c,b ; Ry will also be evaluated for q = t. 

Before presenting the results, a few words have to be said about the parameters entering 
Eqs. (4.2, 4.3). For the masses we will take the values given in Eqs. (3.4, 3.5). For a s (/a 2 ), 
the “natural” scale is fi ~ mC F a s , as shown in Eqs. (3.11). Nevertheless, for c, b we will 
choose n = 2 rn q (but see subsection B for a discussion of this). For tt production, and since 
the choice of /a is less relevant now, we will take /i = Mz ; then we may use directly the 
value of a s deduced from the Z decays. Thus, we take 


a s (M z ) = 0.119 ±0.003 , 


(4.4.a) 


for tt ; and for bb and cc, 


®S ( h ) 


4 tt 


(3 0 In /i 2 /A 2 


1 - 


A In In ji 2 ! A 2 
/3 0 2 In /a 2 /A 2 


(4.4.b) 


Id = 2 m c , 2 rn h ; A = 200+^ MeV , A = 102 - 38 n f /3 


Finally, for ( a s G 2 ) (which does not play a very important role in our evaluations) we choose 
the standard value! 12 ! 
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{a,G 2 ) = 0.042 ± 0.020 GeV 4 . 


( 4 . 5 ) 


As stated several times, our calculations are valid in the nonrelativistic regime, i.e., to 
corrections 0{y 2 ). For the numerical evaluations we will take v < UMax , %ax = 1/2 . 

B. Comparison with experiment: cc and bb. 

The prediction of onr calculation (4.3.a) for cc is shown in Fig. 3. In order to display 
the dependence of the calculation on the choice of the renormalization scale, //, we present 
these RV for two choices of /./: electing /i such that the radiative corrections vanish, i.e., 
such that Co = 0, or fixing // = 2 m c . There is little difference between both choices; here we 
favour this last choice because it ties with the election at higher energies (s 1//2 m c ) where 
one takes fi = s 1//2 . Also shown are the results of a parton model (R y and parton model 
plus order a s (m c ) correction (the last denoted by R x c ^). It is seen that there is a partial 
cancellation of the Fermi factor, 

7T C p Ol s / V 

^ _ g— TrCpOe/v ’ 

and the radiative correction, 2 c 0 (k) (c.f. Eq. (4.3.a)) in such a way that does not 
differ much from f?/ The nonperturbative contribution is only important right at 

threshold. 

Comparison with experiment is shown in Fig. 4. Although the quality of this is not 
enchanting as a fit, a few things must be said in its favour. First, the theoretical curve 
runs, more or less, through the middle of the experimental points. It is clear that the full 
ffV (exact) more centered than the purely partonic f?/ ((y> or what we would have obtained 
neglecting the radiative correction, ( no correc -) (in this last case, the improvement is slight). 
On the average , R% ( exact ) represents a good mean of experiment, which is an interesting fact, 
particularly since part of the dispersion of the experimental points is doubtlessly due to error 
fluctuations. 
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Fig. 5 shows the comparison with experiment for bb. The conclusions are similar to those 
for cc: Rl ^ exact ' ) improves Rl ^ and represents a good average of the experimental points, 
which are now scantier. 


C. tt production. 

We plot in Figs. 6, 7 the quantities R v , R A relevant to production of tt by a vector, axial 
current respectively. It is not possible to use them directly to predict experimental output 
because our formulas do not take into account the width of t. This can be done with the 
standard methods^ 13 !, and we will present the details separately. 

V. THRESHOLD EFFECTS ON “LOW ENERGY” CORRELATORS. 


Consider a correlator, 11(f). It is possible to prove quite generally that it verifies the 
relation (dispersion relation) 

t r°° , Im IBs) 


n(f) - n(o) = - / ds 

IX J 


(5.1) 


s(s — t ) 

By expanding in a power series in a s it follows that Eq. (5.1) is also verified order by order 
in perturbation theory: 

t r°° Im Il( n )(s) 


n (n >(t) -nW(o) = - / ds 

TT J 


(5.2) 


s(s — f) 

When used in this last form, the dispersion relation is little more than a calculational 
device which may simplify the evaluation of hb'd as, generally speaking, Im hb n ) is easier to 
calculate than Re hb n b This method is followed e.g. in Refs. [10] for IIy\ When used 
in the form (5.1), however, a dispersion relation may yield new knowledge. This happens 
when experimental data may be used as input for ImII, thus obtaining values of II in 
regions inaccessible to theory, as is the case for evaluations of the hadronic corrections to 
the anomalous magnetic moment of the muon. Another situation, which is the one we will 
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encounter here, is when nonperturbative methods are employed to evaluate Im II (or parts 
thereof). 

Let us define the threshold effect contributions to 11(f) — 11(0), to be denoted by A^(f), 
to be the quantity 


A th (*) = - / ds 
71 J 

f {n \s) = ImII(s) - 


/ (n) (s) 
s(s — f) 

n 

Y, Imn (n,) (s) 

n '=0 


(5.3) 


That is to say, A^ incorporates the exact contribution to ff(f) —11(0) from a region around 
threshold, up to an energy %, but subtracting the first n terms in perturbation theory. 
Clearly, one has 


n 

nw - n(o) = •£ {n<“'»(«) - n<»'>(o)} + a£>(*) + a£> , (5.4) 

ii'=0 


where A^ would be 


Abb) = 1 fds 

7T J sm 


ImII(s) - Imn (n,) (s) 


n '=0 


s(s — t) 


(5.5) 


Expression (5.4) is then useful when we can argue that A^ is small with respect to A^\ 
or (as happens for fly) when the difference 


n 

Im If(s) - Y Imn (n,) (s) 

n'= o 

may be calculated for s > %. 

It should be clear that Aj:^ + A^ depends on sm; only if exact evaluations were used 
that an exact match (and thus cancellation) of the dependence of each of A^, A^ on s M 
would occur. In favorable cases we expect, however, that the residual dependence would be 
slight. 
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A. The axial correlator. 


We do not have any information on for the axial correlator, but we will give results 

A (2) 

on A ) hA for completeness. Separating the contribution from the bound states and the piece 
above threshold, we write 


A (n) = A yn> 4- A 1 

A ^pole A ' ^a.t. A 


(n) 

+ole A 


(n) 


(5.6) 


Using Eq. (2.4.b) for Ap^ le A , and Eqs. (4.1), (4.2.b), (4.3.b) for A [p [ A , we obtain 


( 2 ) 


A 


(i) 

pole A 


(t) = 


3 N r t 


E 


i 


2 m 2 it ^ Mm Ml h - t 


Kl(0)f 


J N1 ila Nl 

where the sum runs over the t = 1 bound states, and 


(5.7) 


Ai'UW = 


N r t f s M 


ds 


12 7T 2 Js th S (s — t) 
x(u 3 (l-2^ p ) [1 + 2 Cl (k)] (l- C r a * 


4v 2 


nC F a s /v 

_ p—irCpOte/v 


(5.8) 


nC F a s v z 


k = mv , v = J 1 — 4 m 2 j s . 


For the case of ti production, the nonperturbative corrections Sj}\ p are negligible, so in this 
particular case we may approximate the threshold s th by 4 m 2 . Moreover, the spectrum 
of bound states is Coulombic with a very good approximation up to N ~ 4; since bound 
states with N > 4 contribute very little to Eq. (5.7) anyway, we may take the spectrum to 
be Coulombic all the way. Thus, we may rewrite Eqs. (5.7), (5.8) as 
(i) _ N c t ~ 1 N + l 

pole ak 1 ) - ^ E (M^-t) N*(N-l)a(N,l)s ’ ( ‘ - a ) 

a(N, 1) given in Eq. (3.7), and the M N are as in Eq. (3.6) (neglecting the NP piece), 


M N1 = 2 m 11 
Likewise, 

aLYaM = 


C F a s (fi 2 ) 2 C F (3 0 a 2 a s 


8 N 2 


le _ +i r ds 


Sit N 2 


In ^1 + i/}(N + 2) 
rriC F a s 


(5.9.b) 


12 7T 2 J4 m 2 s(s — t) 


X < U 3 [1 + 2 Ci(k)] 1 


C F a s \ nC F a s /v 


"^7-. -V 3 -1 rC F a s v 2 
4v 2 I 1 _ e -nC F ae/ v 


(5.10) 
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The practical interest of this is fairly limited so long as we have no reliable estimate of 
A . We then turn to the vector correlator, for which such an estimate exists. 


B. The vector correlator. 


The calculation is very much like for the axial case, except that we use Eq. (2.4.a) for 
the poles, and Eqs. (4.1), (4.2.a), (4.3.a) above threshold. We End, for tt and neglecting 
the gluon condensate contribution, 


A* 1 ) it\ _ 1 / y'_1_ l 

p ° kv{ ' * h (M& - t) N*a(N, Of ’ 


(5.11) 


and now 


M™ = 


ir N 2 [ mC F a s 


+ ip(N + 1) 


(5.12) 


We split the region above thresholds into two parts, a low energy (he.) and a high energy 
(h.e.) part, according to v < 1/2 or v > 1/2. Note that v = 1/2 occurs for s = 16m 2 /3. 
Later we will discuss the joining of the two regions. 

At low energy we have, 


Hi N t rT6rn 2 /3 

A 1 e ) y (t) = ^4/ ds hcA } , 

12 TT 2 Jim 2 sis — t] 


(5.13.a) 


and, from Eqs. (4.3.a), (4.2.a) 


/£!(»)= Rf - R v t (0)+(I) 

TrC F a s {fi 2 )/v (v(3 — v 2 ) ( 6v 3nv 2 

l __ e -nC F oi3(n 2 )/v | 2 Y TT 4 

v(3 — v 2 ) (3tt 6v ttv 2 \ . 2 . 

2 “ - (l tT H 2~J C ' a A) ’ 

k = rnv , v = J 1 — 4 m 2 js , /j = Mz ■ 


C F a s (ii) \ [1 + 2 c 0 (h)] 


(5.13.b) 


For high energy there exist evaluations^ 14 ! correct to errors (1 — v) 2 a 2 , a/. Subtracting 
from this the order a s piece, we may write the result as 


A 1 ) _ f _ Ao+i) 

J h.e. J h.e. J h.e. 


(5.14.a) 
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where /h° c +1 ^ is the piece of order (0 + 1) in a s 


/d +1> = -f(l-*’) 2 + i(l-0 3 


[9 , 9 / 2 3\ . , 2 1 

[2 V ; 2 V 1 - v 8 ’ \ 


2 


and^ 14 !, 


=-|(i - ++|(i - ®) 3 +[|(i - ”) + \ ( ln 


(1 - v)‘ 


2 


<+( S ) 


3 9 


+ f 3 +-(1 -v)C F 8.7 

V 7T / 2 V 7T 


<+( s ) 


+ 45.3 


Here, 


(5.14.b) 


(5.14.c) 


v—Jl — m 2 (s)/s , r 2 = 1.986 — 0.115 n/ , 


r 3 = -6.637 - 1.2 n f - 0.005 nj - 1.24 ]T Q f 


and the running mass m (/t) is given in terms of m by 


m = m 


a s (+) 

a s (m) 


Aa s (n) - (C F - A)a s (m) 


a _ #7o - A7 i j _ _Jo 

A 2 ’ ™ A ’ 


7o = -3 C F , 7i = — 


3C r 2 97 C f C 4 , 5C F rif 


r A _|_£_ 

6 3 


Finally, 


r * 

12 7+ J16m 2 /3 


2 /3 S(S — t) 


(5.15) 


Ap^jgy, A| 2 y , and A^ y should be compared with the (0 + l)-order direct calculation of 
n: for small t, 




N c t I f 3 t 205 C F a s (n 2 ) 

12 7T 2 5 771 2 1 28 771 2 54 7T 


(5.16) 


Numerically, 
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0.19 ±0.005 


(5.17) 


3 Ml = 1q _ 2 . 205 C F a s {Mp = 

28 m 2 ' ’ 54 vr 

The error in the above equation is due to the experimental error in a s (Mz)- 
As for the A^ 1 / we write, for the various contributions, 

A«(i) = ■ (5-18) 

12 7r z 5 md 

In this way we may compare directly with Eq. (5.17). Then, 

Apoiev(0) = 2.33xl0- 2 

A} 1 ] y (0) = 1.56 x 10 -2 (5.19) 

A« y (0) =2.15 x 10- 2 

The dependence of the A on t is very slight, up to t — M§, where we have 

= 2.52 x 10- 2 

Af.e. v(M|) = 1.67 x 10” 2 (5.20) 

±±(Af!) = 2.24 x 10- 2 

The largest error of the A is due to the error in the mass of the t quark, about which 
little can be done at present. Another source of error is due to extrapolations: we have used 
Eq. (5.13) for /i. e . up to v — 1/2, and Eq. (5.14) for /h. e . down to same value of v. We can 
smooth this rather crude joining of l.e. and h.e. regions by defining 

fl = (1 - v) /l.e. + V /h.e. (5.21.a) 

and integrating this f\ over all the interval, from 4 m 2 to infinity. Eqs. (5.19) are not 
changed substantially; we obtain now, with self-explanatory notation, 

A <1 2+h.«.v',i(0) = 4.1 x 10- 2 . (5.21.b) 

Another possibility is to write 

h = (1 - v 3 ) /l.e. ± V 3 /h.e. , (5.22.a) 
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which respects the terms in v°, v, v 2 which are known at low energy. Then, 

A 1 ( ‘ ) +h ,„.v, 3 (0) = 3.3xl0- 2 . (5.22.b) 

We consider Eqs. (5.21) to be the more reasonable estimate, and take its difference with 
Eqs. (5.19), (5.22) to be a measure of the systematic theoretical errors in our calculation^. 
Thus we finally get 

^inV(O)— ^poL y(0) + Afi y(0) + Ah 1 . v(0) 

= ( 4.1 ±0.8 ) x 10~ 2 . (5.23) 

This may be compared with the evaluations of Refs. [8,9,15], which are clearly improved 
by our results. In particular, the unstabilities noted by Gonzalez-Garcia et alJ 15 l disappear 
almost completely. This is due to our use of information from the high energy region (which 
eliminates the uncertainties due to the dependence on an energy cut-off), and inclusion of 
the radiative corrections, which reduce drastically the arbitrariness of the choice of the scale 
p of a s (/i). 

Although the evaluation is reasonably reliable, it should be clear that the effect is small 
in the sense that A a uy as given by Eq. (5.23) is smaller than the perturbatively known 
piece, Eq. (5.17), by a factor of about 5. 
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APPENDIX A: WAVE FUNCTIONS IN THE CONTINUUM 

We present here some of the technicalities used to solve the radial wave equation in the 
continuum. Let us recall that we are only interested in obtaining the value of the wave 

3 To which one should add errors due to experimental errors in nit, a s (Mz). 
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function at the origin for £ = 0, and its derivative at the origin for £ — 1. This simplifies 
considerably the calculation. The main results of this Appendix were given in Eq. (3.11). 

To End the radial solutions Ru(r) of the Hamiltonian H , as given in Eq. (3.2), we treat 
Hi to first order, and use the method of variation of constants. In other words, we want to 
solve 


Ru \ r ) + ~ Rke ( r ) + 


k = mv , a = 


+ + n u (r) = 0 . 

ar ' 




ar 


A = 


maC F /3 0 a 2 


(A.l) 


mC F a s (p 2 ) ’ An 

taking A infinitesimal. For A = 0, the normalized solution (see Eq. (2.8)) which is regular 
at r = 0 is 


Ru(r) = e iS e^ 2ka e~ lkr (2 kr) 1 M(A e , 2£ + 2, 2 ikr) 

where 


(A.2) 


= £ + 1 + -— , 
ka 

M is the Kummer function, and 5 an arbitrary phase, which we choose to be zero. 
Let us now discuss the case £ = 0. The general solution of Eq. (A.l) is 


(A.3) 


RkO ~ RkO + dRkO i 

5R k o — e K ^ 2ka | T(A 0 ) | e~ ikr Y (r) , 
Y(r) = c 0 (k) M(A 0 , 2, 2 ikr) + 2 ikr 


M(A 0 ,2,2ikr ) / dpipi 
Jo 

+U(A 0 ,2,2ikr) [ dpcp 2 
Jo 


where 


<pi = UX/W , ip 2 = -MX/W , X = ^ M(A 0 , 2, 2 ikr) , 

rCCL ZZrCV 


U is the Kummer function which is singular at r = 0, and W is the Wronskian of U and M, 


W~ 1 = T(A 0 ) e~ 2ikr (2 ikr) 2 
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A term in Y of the form C(k)U is excluded by the condition of regularity at r = 0. 
The constant c o must be determined from the condition of normalization of R k0 , which is 
equivalent here to the orthogonality of R k0 and 5 R k0 : 


dr r 2 R k0 5R k0 = 0 


After matching infinities, the value of co(k) follows: 

T(A*) e w / fco 


Y2 


c 0 (k)— —2ik J drr 2 j <pi — 

= 4< * ArW ( MV + IAAA 

a Jo [ 2 


\Mf 


(A.4) 


The arguments in both Kummer functions M, U are the same, e.g., M = M(Ao,2,2ikr). 
As we stated earlier, we are only interested in i4 0 (0). We then have, 


R k o(0) = e~” /2ka r(A*) [ 1 + co(k) ] 


To complete the calculation it only remains to evaluate Co(k). We will come back to this 
later. 

Repeating the above analysis for the case £ — 1 gives: 

(A.5) 

r =0 

where 


R, 


k 1 


= '1 + Ci{k) } 


R, 


hi 


r=0 


c\(k) = A 


(2 ikf T(Ai 
3 a 


dr r 3 In (fir) M* 


U 


r/ka 


T(AJ) 


12 


M 


(A.6) 


In the above, M, U = M, U(Ai, 4, 2 ikr), and Ai was dehned in Eq. (A.3). We now explain 
in some detail how c 0 (k) and C\(k) are evaluated. 


1. Evaluation of co(k). 


Using the identity: 


M(a, c, 2 ikr) 


r ( c ) c, 2ikr) + £M e ^(«-c) e 2 *r f/ / _ _2i kr ) , 

T(c — a) T(a) 


(A.6) 
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and making a change of variable: r = p/2k, we rewrite Eq. (A.4) as (below r) = ka): 


-i e~ n/v 


Co(fc) — ^ 0 I r / A \ |2 ^0 ) 

2 V | r(A 0 )| 

■/o = J o dr r In (^) e~ ir U(A 0 , 2, i r) 2 r(A 0 ) 2 - c.c. ; (A.7) 

(c.c. = complex conjugate). If we now also rewrite C/ 2 r(A 0 ) 2 as 

U (A 0 , 2 , ir) 2 r(A 0 ) 2 = [f/(A 0 , 2 , ir) 2 r(A 0 ) 2 + ^ , 

then the square bracket can be rotated to r = — ip (p > 0) while its complex conjugate is 
rotated to r = +ip. Eq. (A.7) becomes: 


Jn = 2i / dr r In — 


rp smr 


2 k r 2 


f°° dr r e~ r In ^ \u(A 0 , 2 , r ) 2 r(A 0 ) 2 - c.c.' 
Jo 2k 1 


dr r e 


2, ry r(A 0 ) 2 + U (Aq, 2, ry T(Aq) 2 — 


(A. 8 ) 


Note that each integral is convergent. If we define 


POO 

K 0 (e)= dr r 1+e e~ r U (A 0 , 2, r ) 2 r(A 0 ) 2 , 
Jo 


(A.9) 


(convergent for e > 0 ), we have: 

J 0 = !Ir (lnA- 7E ) - |[/r 0 ( £ )-if o(e )*]^ + 


ln^- \K 0 (e) - K 0 (e)*} + ^ |*r 0 (e) + K 0 (e)* - 2 r(e)l . (A.10) 

2k L J e= 0 + 2 L Je= 0 + 


To evaluate A" 0 (e), we replace one of the C/’s by its integral representation 


e- r C/(a,c, r)r(a) = [' dt e~ r/t t~ c (1 - t ) a " 1 , 
Jo 


(A.ll) 


and use the result 


, 6 -i _ srm v T(6-c+l) _ b ^ r(a + n)r(6 + n) (1 - s 1 ) n , 
/o r(a) r(a + b - c + 1 + n)r(l + n) 


Re(s) > 1/2 . 


(A-12) 


We find: 
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(A.13) 



Jo — ——( 27 e + VKA) + VKA) + (A) — 1 ) 4 1 (A) + (A — 1 ) ^ (A)) 

+ (in ^ - 27^ (mt + A A) - A A) + (A) - 1 ) A(A) - (A - !) E(A)) (A. 16 ) 

00 1 I 77, 

+ E M , „32 ( 2 A A + n + 1 ) - A™ + 2 )) - c.c. . 
n = 0 L + n) 

For small velocities (i.e. for rj = ka < 0.1), we can use the asymptotic behavior of the if} 
function and its derivatives to rewrite Eq. (A. 16) as 

Jo = 2in (inf-1-27, + ^+^ + -..) , (A.17) 

and our result for c 0 (k) is: 

c 0 (k) = - L^- 1 - 27 b + —+ —+ ...1 • (A.18) 

uw 2 [ 2 12 40 

For rj > 0.1, the sum in Eq. (A.16) is evaluated numerically. See Fig. 1 for a plot of c 0 (k) 
as a function of 77 , when /i = 2 /a. 
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2. Evaluation of c\(k) 


For the case l = 1, we need to evaluate 


cAk) = A- tj , 

2h I r(A 1 )| 2 

./, = /°° dr r 3 In (^-) \e~ ir U(A u 4, ir) 2 rfAi) 2 - c.c. 
Jo V2A;/ L 


Ji is evaluated in the same maner as Jq. We first rewrite £/ 2 T(Ai) 2 as 


f/(A 1 ,4,ir) 2 r(Al 1 ) 2 = [c/(A!,4,ir) 2 T^) 2 - V{ir)\ +V(ir) , 
__ 4 12-4 A 2 21 - 16 A l + 3 A 2 




+ 


Then the square bracket is rotated to r = —ip, and its complex conjugate to r 
(A. 19) becomes: 


J 1= Tdrr ln^ 
Jo 2k 

poo 


e ir V{ir) — c.c. 


r/i 


+ J dr re r In — |[/(Ai,4,r) 2 T(Ai) 2 — V(r) — c.c. 

in f°° , _ r 

-/ dr r e 

2 Jo 


[/(A, 4, r) 2 T^) 2 - V(r) + c.c. 


where each integral is convergent. The hrst integral is easily done and gives: 


^ ,, h \ / 3 tt 8 \ . (8 N 


To compute the remaining integrals, we define 


KAe) = / dr r 3+e e“ r U(A U 4, r) 2 T^) 2 , 
Jo 

roo 

LAe) = / dr r 3+e e~ r U(r) , 

Jo 


which are convergent for e > 2. Eq. (16) is then rewritten as: 


Ji = + lim 

e —>0 


(A'lW-A-ior-i.w+iior) 




(A.19) 

(A.20) 
ip. Eq. 

(A.21) 

(A.22) 

(A.23) 

(A.24) 
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Note that the result for e = 0 is obtained by analytical continuation. The results for K\ and 
Li are easily obtained if one uses Eqs. (A. 11) and (A. 12): 


L\(e) = 4 T(-2 + e) + (12 - 4 A{) T(-1 + e) + (21 - 16 A, + 3 A\) T(e) 

r(4 + n + e)r(A! + n) 2 
^ r(l + n) T(Ai + n + 1 + e) 2 


Ad(e) = r(l + e) 2 E 


(A.25) 


To extract the pole part of K\ (e), we rewrite the last equation as: 

T(4 + n + e) T(Ai + n ) 2 


Ai(e) — T(1 + e) 2 E M n(e 

71=0 


e) 2 E 

n =0 


where 


M„(e) 


I' 1 + n) T(Ai + n 
5-2Ai + e(7/2-2Ai) 


- u n (e) 


(A.26) 


(n + l) _1+e (n + l) e 

11 - 12A 1 + 3 Af + e(175/12 - 187R + 5Af) 


+ 


(n + l) 1+l 

The first sum which contains the pole part gives 

'7 

u n (e) — £(—1 + e) + (Re) |5 — 2 4i + e( 

n =0 


(A.27) 


E u n( e ) ~ C ( — 1 + e) + C(e) 
+C (1 + e ) 


5-2A 1 + el--2A 1 


/175 

11 - 12 A! + 3 A 2 + e (-18 A x + 5 A 2 

\ 12 


(A.28) 


In the above equation, ( is the Riemann’s zeta function. Now the second sum in Eq. (A.26) 
converges for e = 0 and so does its derivative with respect to e. We obtain the following 
partial results: 

T(4 + n + e) r(A! + nf 


E 

n =0 L 


T(1 + n) T(Ai + n + 1 + e) 2 


M n (e) 


- M^3 + E ^) + ( 1 + e 



(A.29) 


Q OO 

— E 

cyt n=0 

(2 + 


T(4 + n + e)r(A 1 + ny 


T(1 + n) T(A! + n + 1 + e) 2 
5 

— 2 


- w n (e) - c.c. 


V 


j (V(Ai) - V’(AJ)) - E ^(Ar) +-0(At)^ - E(hi2-- E-8-,/.;) 

^(n + 4) — 2^(Ai + n + 1) 


+* (4 + 1 ) E 

W vj e> 


(Ai + n) 2 


+ c.c. 


(A.30) 


fi + E^ v 

ip(n + 4) — 2il>(Ai + n + 1) 

r* r* 

V W ^0 

C.C. 

Ai + n 
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where we have used A\ = 2 + i/ 77 . For small velocities, we have: 


Ji = 2in 


('“T “ 27e ) ( ! + ?) + ? " 12 + II) 1,2 “ 


1 

7] 7T 


+ 0( V 4 ) , (A.31) 


and our result for Ci(/c), rj < 0 . 1 , is: 

ci{k) = A ^ln ^ — 2'Je 


3 2 4 \ 49 2 167 4 

— h V — V ) — 2 H- ?? 2 -?? 4 

2 V 24 80 1 


+ — (r) - rf + r] 5 ) + 0{rf) 
2 7T v y 


(A.32) 


See Fig. 2 for a plot of c\{k ) as a function of 77 , when 7 / = 2/a. 

Before we close this section, we would like to add a few words about the correctness of 
the results presented here. We have checked that the terms which are proportional to In [i 
in Eqs. (A.19, A.20) are the same as those obtained by solving Eq. (A.l) directly, with 
the obvious replacement In fir —> In fi ; they can be deduced from Eq. (A.2), with r = 0, 
and a -1 —> a -1 + a _1 Aln n . These terms are needed to show that the physical results (see 
for instance Eq. (3.9)) are, to the order we are working, /^-independent (this is of course 
a consequence of the /^-independence of the Hamiltonian H in Eq. (3.2)). We have also 
checked that Eqs. (A.18, A.32) reproduce the correct results in the limit k —» 0. The latter 
are obtained as follows: one first regulates the integrals in Eqs. (A.7, A.19) (e.g., by adding 
a term e~ re in the integrands), then one uses 


lim T(a — c + 1 ) f/(a, c, z/a ) = 2 z 1 ^ 2 c ^ K c -i (2 z l ^ 2 ^j , 


where K c _i is the Bessel function of the second kind. The integrals one obtains involve 
a product of Bessel functions, and are easily done. This provides a crucial test of the 
correctness of the calculation. 
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FIGURES 


FIG. 1. 


Co versus rj = ka. 
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FIG. 2. 


Ci versus r) = ka. 
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FIG. 3. 


v.s. y/s, with m c = 1.57 GeV , A qcd = 0.2 GeV, (a s G 2 ) = 0.042 GeV 4 . The solid 
and dotted lines were obtained with the choice n — 2 m c ; the dashed-line was obtained by 
choosing // such that the radiative corrections vanish (be. such that Co = 0). This shows 
that is not very sensitive to the choice of /j. 
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FIG. 4. 


Theoretical predictions for together with experiment results (Ref. [15], Particle 
Data Book). The dotted-line represents the parton model prediction, and is obtained with 
m c = 1.57 GeV, fi = 2 m c , A qcd = 0.2 GeV, and ( a s G 2 ) = 0.042 GeV 4 . The solid 
line represents the “exact” value of R% and the gray area surrounding it represents the 
uncertainties which are obtained by varying m c , Aqcd and ( a s G 2 ). The dashed-line gives 
R c when radiative and nonperturbative corrections are neglected. The experimental data is 
represented by the diamonds; the horizontal error is less than 4% while the vertical errors 
are large ( rs j 5-35%). 
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FIG. 5. 


RV v.s. y/s. Same conventions as in Fig. 4. 
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FIG. 6. 


Prediction for R) v.s. a/s, with = Mz, a s = 0.119, and ( oi s G 2 ) = 0.042 GeV 4 . 
R] (-exact) j g pitted for m t = 165 GeV (solid line), m t = 180 GeV and m t = 150 GeV 
(dashed lines). 


35 



FIG. 7. 


Prediction for Rf v.s. y/s, with m t = 165 GeV, and (a s G 2 ) = 0.042 GeV 4 . We plot 
Rl (exaci:) using two different criterias for the choice of //: (a) choose fi = Mz and a s = 0.119 
(solid line); (b) choose /i such that the radiative corrections vanish, i.e. such that c\ — 0 
(dashed line). 
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